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Renormalization group equations are derived for the case when both valley splitting and intervalley 
scattering are present in a two- valley system. A third scaling parameter is shown to be relevant 
when the two bands are split but otherwise distinct. The existence of this parameter changes the 
quantitative behavior at finite temperatures, but the qualitative conclusions of the two-parameter 
theory are shown to be unaffected for realistic choice of parameters. 

I. INTRODUCTION 

Renormalization group (RG) studies of multi-valley two-dimensional electron gas (2DEG) systems has been very 
successful in quantitatively describing the transport properties of electrons confined in silicon inversion layers (MOS- 
FETs) ^'12, 13 a disordered medium, for temperatures ksT < h/r, where 1/r is the elastic scattering rate, the 
propagating modes are diffusive, and it is now well understood that these modes play a central role in determining the 
transport properties at low temperatures ^. In two dimensions, in particular, the effects of diffusion are profound. The 
electron-electron (e-e) scattering amplitudes, for example, develop non-analytic corrections that result in enhanced 
correlations at low energies ^ . It has been shown that RG theory applied to a weakly disordered system is able to 
capture this scale (energy or temperature) dependence to all orders in the e-e scattering amplitudes, making it the 
most promising analytical technique available to understand the physics of disordered systems. (Pedagogical reviews 
of the RG theory can be found in Ref. ^'^.) 

Weak disorder implies that h/r < Ep-, where Ep is the Fermi-energy. Typical high mobility two-dimensional 
semiconducting devices have very small Fermi-energies with a scattering rate which is even smaller due to the very 
high mobility of the samples making it very difficult to access the diffusive region at experimentally reasonable 
temperatures. Si-MOSFETs on the other hand have only moderately high mobilities so that H/t is of the order of 
a kelvin while Ep is of the order of a few kelvin. The impurity scattering in these inversion layers is short-ranged 
in character making quantum scattering the dominant scattering mechanism, while semi-classical effects arising from 
the impurity potential landscape are negligible at low temperatures. For these reasons, as noted in the beginning, RG 
theory has been particularly successful in describing the properties of electrons in silicon inversion layers. (See Ref. ^ 
for how the diffusive regime is identified experimentally and for a quantitative comparison of theory with experiment.) 

The conduction band of an n-(OOl) silicon inversion layer has two almost degenerate valleys located close to the 
X-points in the Brillouin zone ^. The abrupt change in the potential at the interface, which breaks the symmetry in 
the z-direction perpendicular to the 2D plane, leads to the splitting of the two valleys. Although intervalley scattering 
originates from both impurity scattering and scattering due to e-e interactions, the imperfections at the interface, 
which are distributed on the atomic scale, are the main source of the large momentum transfer Qo in the z-direction 
needed for intervalley scattering. 

The RG theory developed in Ref. considered the valley degrees of freedom to be degenerate and distinct, hence 
quantitative comparisons with experiments performed in Ref. ^ were limited to temperatures larger than the valley 
splitting, A^, and the intervalley scattering rate, h/r^^ both scales being sample dependent. This paper develops the 
relevant scaling equations in the presence of valley splitting and intervalley scattering. 

The scaling equations are presented in three different temperature regimes: (i) high temperature region, T > Ty 
and T*, where ksTy = Ay and ksT:^ = ^/^x, (ii) low temperature region, T <T^, (iii) and intermediate temperature 
region, <T<Ty. The last of the three regions is relevant when the band splitting is large so that effective mixing 
of the valleys due to impurity scattering occurs only at sufficiently low temperatures; it is shown that the standard 
two-parameter description has to be modified in this case to include a third scaling variable which has quantitative 
effects at finite temperature but does not affect the asymptotic conclusions of the two-parameter theory. 

II. DIFFUSION MODES AND FERMI-LIQUID AMPLITUDES 

Electrons in valleys can be conveniently labeled using additional valley indices = ±. (For our purpose, the 
number of valleys Uy = 2 located at iQo^? where Qo ^ 0.85 x (27r/a) with a being the lattice constant of silicon.) 
This increases the number of single particle states to (spin) x (valley) =4. Since the diffusion modes, responsible for 
the relaxation of density and spin perturbations (and valley in our case) in a disordered system at long times and 
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distances, formally occur via particle-hole excitations, the corresponding number of (particle) (g) (hole) diffusion modes 
equals 16. This is a four fold increase from the case of one valley and has significant quantitative effects on transport 
as shown in Ref. At low temperatures, some of these modes develop gaps (cut-offs) proportional to and A^^, 
and are therefore ineffective (non-singular) for T below the characteristic temperature scales Ty and 10,11^ leading 
to quantitatively different scaling as the temperature is varied. 



A. Single particle properties 



FIG. 1: Diagrams contributing to l/r in Eq. (2) are shown. The properties of the intra- and intervalley impurity scatterings 
are defined in Eq. (1), represented here by single and double dashed lines, respectively. 



At low electron densities the mobility of a 2DEG is determined by the charged centers within the Si02 layer. Due 
to the short-ranged nature of the impurity scattering in silicon inversion layer structures, the Drude relation for the 
mobility, ji = er jm^ gives a direct measure of the single particle life-time, r. Here, e and m are the charge and 
the effective mass of the electron, respectively. Ando ^ argued that the mobility is also determined partially by the 
intervalley scattering rate. To this end, the two different scattering rates, that is, the intravalley and intervalley rates, 
can be incorporated by introducing two scattering potentials i^(q) and '^(q), respectively. The potential it(q) is 
slowly varying on the scale of 1/a if the impurities in the oxide layer is uniformly distributed, while ^'(q) is a rapidly 
oscillating function with momentum of the order of 1/a. Hence the random average of the potentials ('?i(q)v(q)) = 0, 
with ii(q) and 'L'(q) satisfying 

(^(q)^(qO) = ^q+q'^:^^ (la) 

(^(q)^(q)) = (lb) 
where = m/27r is the density of states per spin and valley. The total life time, r, then equals (see Fig. 1) 

111 

- = - + — (2) 

r rii r_L 



B. Particle-hole diffusion propagators 



The form of the particle- hole propagators (diffusons) for the impurity model defined in Eq. (1) have been calculated 
in Ref. ^^■>^^ , The calculations are extended here to include valley splitting. 

The fiuctuations in the diffuson channel, P((7, uj\ have a diffusive singularity I^(<7, uj) = l/^D^q^ + 1^1)- Finite valley 
splitting and intervalley scattering introduces gaps in V{q^uo) thus cutting off the singularity. The different diffuson 
modes involving fiuctuations in the valley occupations are shown in Fig. 2. The details of their derivation are given 
in Appendix A. 




FIG. 2: The various diffuson blocks in the presence of valley splitting and intervalley scattering are shown. 
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We start by defining tiie elementary diffuson blocks, V\\^u-,^\\,v and 2^±,t^ shown in Fig. 2. The diffuson blocks 
V\\ are insensitive to valley splitting since both the particle and the hole (corresponding to the top and bottom lines 
with arrows moving to the right and left, respectively) belong to the same valley. The index for the ^^±,t^ diffuson 
indicates the valley index of the particle, with the hole being in the — valley; the two valleys are nonequivalent for 
finite A^. In Appendix A, the equations satisfied by the diffuson propagators are solved in the limit of weak splitting 
A^r < 1. The solutions are expressed in terms of the diffusons V± = V\\^u i and ^^±,t^, with the corresponding 
gaps A_ = 2r||/r(rx — t\\) and A^ = t\\/tt^. (Note that corresponding to the valley "singlet" mode, is gapless, 
hence A+ = 0.) 

In the limit when the intervalley scattering is much weaker than the intravalley scattering, i.e., t±_ ^ ry, the 
scattering time r ^ t\\. The gaps in this limit correspond to A_ ^ 2/t±_ and Ax ^ l/r^. In this weak scattering 
limit, relevant to high-mobility MOSFETs, the form of the diffusons obtained in Eqs. (A3) and (A4) reduce to: (the 
overall factor 1/2t:vt'^ is suppressed) 



1 
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L)o^^ + |c^|+2A* 
1 

Doq^ + \u\ -ZT^A^ + A* 



(3a) 
(3b) 
(3c) 



where A* = l/r^. The number of modes that are effectively gapless depends on the relative magnitude of T (or 
frequency) with respect to the corresponding temperature scales Ty and T*. At high-T all modes are gapless, while 
at the lowest T only P+ remains gapless. 



C. Electron-electron scattering amplitudes 
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FIG. 3: The scattering amplitudes including the valley degrees of freedom are classified in terms of the standard static Fermi- 
liquid amplitudes Fi and F2. The same subscript convention used to classify the diffusons in Fig. 2 are used here. 



In this section the relevant e-e interaction scattering amplitudes are identified. These amplitudes are conventionally 
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described by the standard static Fermi-liquid amplitudes Fi and F2 defined in terms of the spin texture of the scattering 
of the particle- hole pairs. The amplitudes are easily generalized to include the valley degrees of freedom. They are 
shown in Fig. 3. Note that the intervalley scattering amplitudes Fix and ^2,\\,v generally negligibly small in a 
clean system because the Coulomb scattering involving large momentum Qo in the z-direction is suppressed when 
the width of the inversion layer is many times larger than the lattice spacing. It is more convenient to work in the 
same basis as that used for the diffusons, i.e., Fi± = | (ri,||,n y and F2± = (r2,||,n ^^2,\\,v)^ as it allows for 
the amplitudes to be easily combined with the diffusion modes. 

III. DIFFUSION CORRECTIONS 

It is now well understood that while the diffusion propagators when combined with e-e scattering lead to the 
appearance of logarithmic corrections to the resistivity (Altshuler-Aronov corrections), the e-e scattering amplitudes 
themselves develop logarithmic corrections due to the slow diffusive relaxation ^. In this section, these logarithmic 
corrections are obtained self-consistently in the limit of weak valley splitting (A^r < 1) and weak intervalley scattering 
(rx > ry). 

T.<r- r, — — ' ~ T. 



FIG. 4: The skeleton diagrams for S are listed. F is the interaction matrix and the hashed blocks are diffusion propagators. 
The 11+ matrix is obtained by adding the contributions from = iT^. The results are presented in Eq. (Bl). 

The e-e interaction corrections to the diffusion propagators are expressed in terms of the "self-energy" matrix S. 
The relevant diagrams are shown in Fig. 4. Expanding to order and uo one obtains, for example, for the 

gapless P+ propagator, the renormalized propagator V^^{q^uj) = Dq^ + zuj^ where D is the renormalized diffusion 
constant and z is the frequency renormalization parameter that determines the change in the relative scaling of the 
frequency with respect to the length scale ^'^ {z = 1 for non-interacting electrons). The corresponding corrections to 
D and z obtained by evaluating the diagrams in Fig. 4 are given in Eq. (Bl) in Appendix B. 

The skeleton-diagrams representing the diffusion corrections to the e-e scattering amplitudes are shown in Fig. 5. 




FIG. 5: Skeleton diagrams for 5Ti^a- By appropriately choosing the F vertices for given values of t'z — =bT^ all the corrections 
^F^,cK, where i = 1,2 and a = ±, are calculated. Each combination of the valley indices comes together with the appropriate 
diffuson matrix elements. Note that going down the rows increases the number of F vertices. 
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(For a detailed discussion of these corrections, see Refs. ^'^.) The calculations are generalized here to include valleys. 
By appropriately choosing the F vertices for given values of = ix^ in Fig. 5 all the corrections, (^F^^c^, to the 
scattering amplitudes F^ q,, where i = 1,2 and a = ±, can be calculated. For example, to calculate (^F^+, since 
Fi+ = F^ 11^^ + ^i,\\,v', the contributions from = iT^ are added, while they are subtracted when calculating (5F^_. 
The results are given in Eq. (B2) in Appendix B. (The corrections to the amplitude ^ are not given as they are 
equal to (5F2+ for T >Ty and irrelevant for T <Ty due to the gap.) 

The corrections SD^ Sz and STi^^ in Eqs. (Bl) and (B2) include all modes, both gapped and gapless. Clearly, only 
modes that are effectively gapless lead to logarithmically divergent corrections. Since the frequency integrations range 
from T < uj < 1/r (the upper cut-off follows from taking the diffusion limit), for T <T^^ both V_ and Vj_ are gapped, 
while only the P_ modes are effectively gapless when < T < Ty. (The P+ mode is always gapless.) Of course, 
when T > Ty and T^, all modes are gapless. As as result, the corrections are clearly sensitive to the temperature 
range considered. 

A. High temperature range: T > Ty and T* 

For T >Ty and T*, all the modes {a = ±, 1.) appearing in Eqs. (Bl) and (B2) are effectively gapless, i.e., they 
take the form V{q^uo) = l/{Dq^ + zoo). As noted below, not all amplitudes Vi^oc are relevant at these temperatures. 
For instance, since intervalley scattering is irrelevant for T >T^^ the amplitudes Fix and ^2,\\,v^ whose initial values 
are vanishingly small, can be set to zero. As a result (see Fig. 5), Fi^ ^ and F2+ ~ F2-. Further more, since valley 
splitting can be ignored for T >Ty, the amplitudes Fi and F2± are indistinguishable from the amplitudes ^i^\\^v 
and F2+, respectively, implying that the initial value of Fi_ = and F2± = F2+. 

It can be seen from Eq. (B2) that choosing the above initial conditions, namely, Fi^ = Fi_ = 0, and setting all 
the amplitudes to be equal, and all the propagators to be gapless, gives ^Fi_ = and ST2- = (^F2+, which 
are consistent with the choice of the initial conditions. Hence, Eqs. (Bl) and (B2) reduce to the form (with the 
substitution F2Q; = F2 and Va = V): 

f = -l JI'^{r,^-^r,)vH,,.)Dq^ (4a) 
= -^/(0(r.,-4r.)p(,,o) (4b) 
(5ri+ = — / ^r2P(<z,o) + 4*(r2) (4c) 

ST2 = — [ -^ri+P(g,0) + 16M/(r2) (4d) 
where // = / dq^/{27T)'^ J du;/{27r) and ^(F2) equals 

The above equations were first obtained in Ref. they correspond to the case when the two valleys are distinct 
and degenerate. 

B. Low temperature range: T < T* 

When T < T*, both V_ and Vj_ are gapped and therefore irrelevant. Hence, only the P+ mode survives. Dropping 
the contributions of the gapped modes in Eqs. (Bl) and (B2) lead to a self contained set of equations involving only 
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the amplitudes STi+ and ST2+- The equations, after dropping the + sign in and r2+, reduce to 



SD 

~d' 


4 1 




(6a) 


Sz 


1 


/ (27rV r2)P(q,0) 


(6b) 




Attu J 




(6c) 


= 


-/ 

TTU J 


'^ri+p(q,o) + 4*(r2) 


(6d) 



These equations correspond to the case when the two vaheys appear as a single vahey due to intervalley scattering. 
(Note that valley splitting is irrelevant in this case as the propagator is always gapped when T <T^, irrespective 
ofT,.) 



C. Intermediate temperature range: < T < Ty 

This limit when the valley splitting is large, so that the intervalley scattering rate <C T^, is interesting. For 
temperatures in the intermediate range < T < T^, only the V_\_ mode is gapped, while both V± are gapless. 
Although the initial value of ri_ ^ when T > Ty (see discussion in Sec. Ill A) it can be seen from Eq. (B2) that 
^Fi- 7^ when T < Ty and is therefore generated at intermediate temperatures. This introduces a third relevant 
scaling parameter distinct from the high and low temperature regimes. (Since T > T*, F2+ = F2-, but because 
T <Ty the F2± amplitude is irrelevant.) 

Dropping the V± terms in Eqs. (Bl) and (B2) and setting F2+ = F2- = F2 and V± = V gives 

-^JJ {Ti.+Ti+-2T2)V^q,u;)Dq^ (7a) 
-;^/(0(r-+r--^^^)^('^'O) 

\ 1 7fT2(ri+ + ri-)P(9,o) + 8*(r2) (7d) 

Note that although both (^Fi+ and 5Ti_ are equal, their initial values are different. 

The relevance of the Fi_ amplitude in the temperature range < T < is specific to problems with split-bands, 
and was first discussed in Ref. ^ for the case of spin-splitting in a multi-valley system. 



5D _ 

Sz = 

sr2 = 



IV. RENORMALIZATION GROUP EQUATIONS 

In Sec. Ill, the leading logarithmic corrections in all the different temperature ranges have been listed. It is now 
possible to set up the scaling equations. To this end, first note that all the corrections involve only one momentum 
integration, and since every momentum integration generates a factor of which by Einstein's relation is propor- 
tional to the resistance p, the corrections are limited to the first order in resistance (disorder). The limitation on 
the number of momentum integrations also constraints the number of e-e vertices in the skeleton diagrams shown 
in Figs. 4 and 5. These corrections can now be extended to all orders in F (but still first order in p) by performing 
ladder summations as shown in Fig. 6. It amounts to replacing the static amplitudes F by the dynamical amplitudes 
U (g, uj) as discussed below. 

Since the ladder summations do not introduce additional momentum integrations, the resummation allows the 
corrections to be evaluated to infinite order in the interaction amplitude leaving p as the only expansion parameter 
in the theory ^ . 

For the amplitudes F2q;, the ladder sums are most easily done in the basis a = ± and ±, as it can be checked by 
inspection that the indices are conserved in the ladder. Using F2q; and Va in Fig. 6, one obtains the corresponding 
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FIG. 6: Extending the static amplitudes F by the corresponding dynamic amphtudes U through ladder summations. 



dynamical amplitude U2a{Q^^)^ where 



The propagators Va are defined in Eq. (3) and 



Dq^^{z^r2a)oo^^c 



(8) 



(9) 



It should be noted in, for example, Fig. 5, that only those interaction vertices involving frequency integrations 

can be extended to include dynamical effects. For convenience, the corresponding F2 vertices are enclosed in square 

brackets in the function ^ in Eq. (5). Substituting for F2 in Eq. (5) with U2 from Eq. (8) (the a index is dropped 

since only gapless modes have been retained in (5)) and performing the q and co integrals leads to the very simple 
6,7. 



expression 



vi>(r2)=('^)x|iog 



1 

TV 



(10) 



The dimensional resistance p = l/4:{27r'^uD) corresponds to (e^/7r/i)i?n, where Ru is the sheet resistance. The factor 
4 arises due to the spin and valley degrees of freedom and u is the density of states per spin and valley. Also note 
that up to logarithmic accuracy the upper cut-off can be replaced with 1/r. Since, the remaining integrals in Eqs. (4) 
to (7) are of the form J dPqV{q^{))^ they can be evaluated directly as 



(fq 



,PfeO)=2plog 



1 

TV 



(11) 



It remains to evaluate the integrals for 3D and 3z. The 3z integrals do not involve frequency integrations and can 
therefore be evaluated using Eq. (11). The 5D corrections, however, contain frequency integrals, and therefore the 
Vi± amplitudes, in addition to F2, are also to be extended to all orders via the ladder sum. 

This is most easily done in the spin-singlet basis 



r.± = r 



i± 



- 2± 



(12) 



This is so, because the spin and valley of the electron-hole pairs in the singlet and triplet basis are individually 
conserved in the ladder sum. (Note that the amplitude is written in the (spin-singlet)(8)( valley- singlet) basis, while 
the ' — ' amplitude is in the (spin-singlet) (valley-triplet) basis; the valley-triplet corresponds to = 1, 6'^ = 0).) The 
corresponding dynamical amplitudes Us±{q-,oo) on performing the ladder sum gives 



Us±{q,uj) = Ts± 



Vs±{q,^) 
V±{q,u) 



where 



(13) 



2^s±(9,w) 



L>g2 + (^-4r,±)a; + A± 



(14) 



(Note that is introduced for notational uniformity, in fact = 0.) 

Special attention is to be paid to the ladder sums involving F5+ when Coulomb interactions are present. The 
Fi+ amplitudes in this case includes amplitudes of the kind shown in Fig. 7, which can be separated by cutting the 
statically screened long-ranged Coulomb line once. They are denoted here as ro+. This distinction is important 
because the polarization operator, I\.{q^uo)^ which is irreducible to cutting a Coulomb line does not include ro+. (The 




FIG. 7: The statically screened long range part of the Coulomb interaction, which can be separated by cutting just one Coulomb 
line is shown. The amplitude ro+ is obtained by adding the amplitudes with t'^ = =b. The shaded triangles represent the static 
vertex corrections V. 

corresponding Fq- and To± amplitudes are zero. The former is identically zero, while the latter involving intervalley 
scattering is vanishingly small.) 

Analyzing the polarization operator, Il{q,uj), provides key insights into the relationship between the various ampli- 
tudes (and z) The form of n(g,a;) is analyzed here in the presence of valleys. In the limit of q^co 0, it can be 
shown that Il{q^u;) takes the form: 

It is important to note that only the F^^ amplitude, corresponding to the singlet mode, appears in the expression for 
Il{q, 00). The factor dn/dji is the thermodynamic density of states and the parameter V is the static vertex corrections 
represented as shaded triangles in Fig. 7. 

The two terms in Eq. (15) correspond to the static and the dynamical contributions, respectively. By construction, 
the static limit Ii{q ^ 0,a; = 0) = —dn/dji is satisfied. In the opposite limit, local conservation law requires that 
Yl{q = 0,cj ^ 0) = 0. From Eq. (15) it can be seen that for the latter condition to be satisfied the following relation 
must hold: 

8n ^^^^ 



dii z- 4F, 



in which case, Ii{q^ uo) takes the form: 



n(^,c^) = X ^ , ^ , (17) 



When Eq. (16) is combined with the definition of Fo+ as the static limit of the Coulomb interaction, i.e., Fo+ = 
V'^dfi/dn, the following expression for Fo+ is obtained: Fo+ = \{z — 4Fs). Hence, conservation laws provide the very 
important relation 

^ = 4(ro+ + ri+)-r2+^4rff (is) 

where F^_^ = Ff:^ — F2+/4, denotes the singlet amplitude in the presence of long-ranged Coulomb interactions. Since, 
only the Coulomb case is considered in the following, all the Fi+ amplitudes appearing in Eqs. (4) to (7) are to be 
replaced by their long-ranged counterparts 

ri+ rf^ = ro+ + ri+ (19) 

Direct inspection of Eqs. (4) to (7) shows that the singlet combination in Eq. (18) is satisfied everywhere, i.e., 
6{z — 4Tg^) = 0, provided (5Fo+ = 0. (This is a well established result, with great importance for the general 
structure of the theory ^'^.) In particular, the corresponding dynamical amplitude U^_^{q^u) reads 

Note that unlike the Us-\- amplitude (see 13) is a universal amplitude independent of r^_^. This is a direct 

consequence of the singlet relation (18) . 

The scaling equations discussed below are obtained from Eqs. (4) to (7) after (i) rearranging all the Ti± ampli- 
tudes to give Fgi, and then replacing F5+ with r^_^, (ii) replacing the static amplitudes where applicable by the 
corresponding dynamical amplitudes and (iii) substituting p = 1/ A{2t:'^vD). 
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It is convenient to express the equation for p in terms of the scahng variables, 72 = V2/Z and 7^ = —AVs-/z. In 
terms of these variables, the equations for p, 72 and 7^ form a closed set of equations independent of z. The final 
RG equations, along with the equations for z^ are given below. The scale ^ = log(l/Tr) is used in these equations. 
To logarithmic accuracy 1/r can be used as the upper cut-off. The range of applicability of ^ is defined in each case 
separately. 



1. High temperature limit: T >Ty and 



dp 



2 

d\nz p 



(1 - 15*(72)) (21a) 
(1 + 72)' (21b) 



(1 - 1572) (21c) 



2. Low temperature limit: T <T^ 



dp 



2 

d\nz p 
di " ~2 

3. Intermediate temperature limit: T^, <T <Ty 



p^ (1 - 3*(72)) (22a) 
(I+72)' (22b) 



(1 - 372) (22c) 



dp 



p' (1 - $(7„) - 6$(72)) (23a) 



^ = f [(l+72)' + (l+72)(72 - 7«)] (23b) 
^ = |(l+7„)(l-7„- 672) (23c) 
^ = -f (1-7. -672) (23d) 



The variable $(7) is defined as 

$(7)= (^l + ^^log(l+7)-l (24) 

The factors 15 and 3 appearing in Eqs. (21) and (22), respectively, correspond to the number of effective triplet 
modes. In the case of two distinct, degenerate valleys, the 16 spin-valley modes break up into 1 singlet and 15 "triplet" 
modes, while in the limit of strong intervalley scattering the two valleys are effectively combined into a single valley 
leading to 3 spin-triplet modes. 

When the valleys are split, as in Eq. (23), the amplitude 7^ plays a significant role as the temperature is reduced 
well below Ty. Given that 7^ = (r2 — 4ri_)/z and that ri_ for T > T^, it follows that 7^ 72 as T approaches 
Ty from above. When T <^ Ty^ the two amplitudes 72 and 7^ diverge from each other significantly. For T < Ty, 
however, it is reasonable to assume that 7^ ^72- This is relevant if the lower cut-off is not much smaller than 
Ty. In this case, the equation for p and 72 pertaining to the different temperature ranges can be combined to give 
dp/d^ = p^(l — {4K — 1)^(72)), and d^2/d^ = p(l + 72)^/2. Here, K = nl = 4. when the valleys are degenerate and 
distinct (high temperature), K = = 1 when intervalley scattering is strong (low temperature) and K = Uy = 2 
when the valleys are distinct but split so that each valley contributes independently (intermediate temperature). For 
direct comparison with experiments, these simplified equations should suffice for most samples. 

The situation changes, however, once T <C T^, but still greater than T^. We see that 7^ and 72 evolve differently 
until jy reaches the fixed point value of 7* = —1, at which point <^>(— 1) = —1. (This fixed point is relevant only 
when T:^ ^ 0.) The system at this point reduces to a single valley system with resistance 2p. The above properties 
are generic to systems with split bands (spin and valley) as has been discussed in detail in Ref. ^. 
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To summarize, RG equations have been obtained in the case when both vahey sphtting and intervahey scattering 
are present. The results can be directly used to compare with experiments in a two- valley system after adding the 
weak-localization contributions, which are not included here. The case when the two bands are split but otherwise 
distinct is quantitatively different due to the existence of a third relevant scaling parameter. The asymptotic metallic 
behavior is, however, not affected. 
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APPENDIX A: DIFFUSION PROPAGATORS 




FIG. 8: Ladder diagrams for the various diffuson blocks are shown above. As in Fig. 1, the intra and intervalley impurity 
scatterings are denoted by single and double dashed lines, respectively. The particle and hole Greens functions for finite 
correspond to Qr^ (/c, e) = [ie — {^k — T2A^/2) -\- i/2T sgne]"-^ 



The ladder diagrams for each of the diffuson blocks, Py ^^, Py ^^ and ^^±,t^5 are detailed in Fig. 8. The corresponding 
equations are given in Eq. (Al). Note that the diffusons are coupled in the presence of intervalley scattering. For 
convenience, the scattering rates in Eq. (1) are defined as Ay = l/27ri/ry and A_\_ = l/27rur±. 

Dy,, = Ay +AyXyPy,, + A^XyDy,, (Ala) 
Dy,, = A^+AyXyPy,, + A^XyDy,, (Alb) 
^±,T. = Ay + AyX^,^^Dx,T. (Ale) 

where 

X\\{q,ij) = ^ (k + q, e + uj)Qt, (k, e) 

k 

^ 2iivT [l-Tuo- rDoq^ + • • • ] (A2a) 
X^^^^{q,uj) = ^ (k + q, 6 + Cc;)g_^^ (k, e) 

k 

^ 27rur [l-rco- tDoq^ + zT^(A^r)] (A2b) 

Here D is the diffusion constant. In the diffusion approximation, i.e., for e(e + a;) < 0, it is sufficient to evaluate X 
in the long wavelength and small frequency limit. Only the weak splitting A^r < 1 limit is considered here. 
Eqs. (Ala) and (Alb) are easily decoupled by defining 

P± = Dy,, ± Py,, = X 2 / (^3) 

where A+ = and A_ = 2r\\/r{r± — Ty). Note that P+ is gapless. Substituting (A2b) in (Ale) gives for T>j_^r^^ 

^^'^^ ~ 2^n^ ^ Doq'^ ^00 -ir^Ay^ ^^^^ 



where A^ = r\\/rr±. 



11 

APPENDIX B: DIFFUSION CORRECTIONS 

The expressions for 5D and 5z extracted from the diagrams in Fig. 4 for the P+ propagator are given below. 



5D _ 4 



Dq^ [Dliq^co) {T,^ - r2+) + Vl{q,u;) {T,_ - T2-) + Vl{q,u;) {T,^ - 2T2±)] (Bla) 
TTU J (27r)2 



Sz = — / --4^[p+(g,o)(ri+-r2+) + p-feo)(ri_-r2-) + Pxfeo)(ri^-2r2x)] (Bib) 



(Note the factor of two in front of T2±'-, for convenience the index is suppressed in the ± terms.) The diffusion 
corrections to the amphtudes F^^c^, where i = 1, 2 and a = ±, are detailed below. 

= i / (0[r2.i>. + r._p_ + 2r.,pj 

^ ^ ^ L0T2+ [Tl+Vl + Tl_V^_ + 2Tl^Vl] - [T2+VI + T2-VI. + 2r2±Pl] 

c.2r2+ [r^^Pj + Tl_Vt + 2r2^Pi] (B2a) 



1 

"2i^ 



^r2+ = ^ / 7^[ri+D+ + Fi_P_+Fi^Dj 



t:v J (27r)2 

r2+ [r2+Pi + r2-p^ + 2r2±Pi] - \ [rl+vl + tI_vI + 2tIj)1] 

4 
2 



^ I j uj^Tl^ [Tl^VX + Tl_Vt + 2rLl>l] (B2b) 

^r,_ = ^/|^[r2,p_ + r2_p,-2r2.pj 

+^ r2- [(r2+ + r2-)p+p- - 2r2iPi] - \ [T2-T2+ {vl + vl) - 2tIj)1] 

ri_ [(ri+ - r2+) (p+p_ - vl) + (ri_ - r2_) (p+p_ -vl)-2 (Fi^ - 2r2i) vl] 
-\ jj wr2- [r2_r2+ (p?.p+ + v_vl) - 2r2_r2±Pl] - ujtI_ [r2+p^p_ + T2-V+v'i - 2r2_Lpi] 
jj uj^rl_ [r2_r2+p^p^ - rl^vl] (B2c) 
^r2_ = ;^/(0[ri.i^- + r._p,-r,,pj 

+^JJ r2_ [r2+p2 + r2_p2 + 2r2iPi] - ^ [r2_r2+p^ + Tl_vi + 2rl^vl] 
11^2- [ri+ (p+p_ - vl) + ri_ (p+p_ - v^_) - 2Ti^vl] 

JT2- [rl+vl + Tl_vl + 2r^^Pi] - wr2_r2+ [r2+p;t + r2_p^ + 2r2±Pi] 

iu'^T2-T2+ [Tl+VX + Tl_Vt + 2ri^Pl] (B2d) 



The terms above are ordered in correspondence with the diagrams appearing in Fig. 5. The square brackets 
gather vertices that come together with the diffuson propagators. Note that the first term is unique in that it 
does not involve frequency integration. (If these equations are calculated using perturbation theory, an additional 
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wave-function renormalization term ( appears ^. The renormalized amplitudes given below correspond to F^^, with 
SC = -I JJ (ri+ - T2+)Vl + (ri_ - r2-)P?_ + (ri^ - 2T2±)V^^. it should be noted that the term C does not 
appear in the non-linear sigma model approach developed in Ref. ''.) 
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